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NOTATION 

Cx,  Cy,  Ci  aero-normal  coefficient  of  the  force  X,  Y or  Z given  by 

Cx  = Xj(\p'Ve*S\  for  example 

Cm  aero-normal  coefficient  of  the  moment  m 

= m/apeVSSIo) 

Cx*.  etc.  derivatives  based  on  aero-normal  force  coefficients  and  aero-normal  variables 

_ 3C,  _ 1 DX  1 DX 

5a  _ beV.*S  is  ~ lpeVeloSU 

Cm±,  etc.  derivatives  based  on  aero-normal  moment  coefficients  and  aero-normal 
variables 

_ 5Cm  ^ 1 1 Jm 

_ 1*"  “ Ipeivsio  J*  = Ip'VeSlo Si 
Iy  moment  of  inertia  about  a y-axis 

j V-1 

K a constant,  derivatives  used  as  spring  rates,  K„  & Kg 

lo  representative  length 

m pitching  moment  about  a y-axis 

mi  mass  in  stability  axes 

trie  = p aero-normal  mass 

mx  pitching  moment  incidence  derivative 

q pitch  rate;  equal  to  6 

S representative  area 

T Tension  or  Torque;  amplitude  is  T0 

V,  datum  velocity 

“ incidence  angle,  usually  sine  definition  as  per  Figure  1 

A change 

r)  pitch  motivator  deflection 

8 pitch  angle 

8C  cable  pitch  angle 

M density  parameter,  aero-normal  mass 

A real  part  of  an  exponential  solution 

pe  datum  air  density 

<u  frequency  in  radians/second 

*°o  natural  frequency,  t»  for  undamped  system 

Superscript  bar  means  aero-normal  terms 
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1.  INTRODUCTION 


To  measure  the  forces  and  moments  on  a manoeuvring  aircraft  by  taking  scale  model 
measurements  in  a wind-tunnel,  it  is  necessary  to  know  the  relationship  between  the  equations 
of  motion  for  the  free  flight  aircraft  and  the  wind-tunnel  model.  The  reasons  for  such  measure- 
ments were  discussed'8*1  in  the  first  Note  of  this  Series. 

The  wind-tunnel  model  system  normally  will  have  additional  restraints  and  constraints  on 
its  flight.  This  note  will  investigate  the  equations  of  motion  and  resulting  forces  for  several 
such  restraints  starting  with  a simplified  version  of  the  free-flight  arrangement.  It  is  hoped 
that  this  will  shed  light  on  the  advantages  and  disadvantages  of  the  various  methods  of  oscillating 
test  models  which  will  be  studied  in  the  next  Note  of  this  Series. 

This  Note  concludes  with  a bibliography  of  introductory  texts.  By  no  means  an  exhaustive 
list,  this  covers  publications  before  1970  which  are  considered  complementary  to  this  series  of 
Notes.  Lists  of  later  publications  can  be  easily  obtained  from  any  of  the  computer  abstracting 
services. 

A single  degree  of  freedom  system  will  be  considered  throughout  and  for  convenience  it 
will  always  be  a pitch  system.  A yaw  system  would  be  identical,  with  the  appropriate  different 
moment,  force  and  inertia  terms,  apart  from  the  lack  of  weight.  In  fact,  any  system  which  measures 
pitch  derivatives  can  measure  yaw  derivatives  (by,  for  example,  rolling  the  model  ninety  degrees). 
There  may  be  difficulties  with  freely  suspended  and  cable  towed  models  if  all  possible  modes 
are  not  stable,  but  in  general  a sensible  choice  of  pivot  location  will  solve  these  problems. 

1.1  The  Aero-Normal  Notation 

For  convenience  and  simplicity  of  equations  and  relationships  between  quantities,  most 
aerodynamic  results  are  expressed  in  non-dimensional  terms.  Several  chemes  for  making 
quantities  “normal”  exist  and  the  ones  of  major  interest  were  consk-red  in  detail  in  an  earlier 
publication  in  this  series. <8*’  The  system  used  throughout  this  publication  is  called  the  aero- 
normal  notation  and  is  identical  to  the  usual  “C”  notation  for  a static  case.  Coefficients  are 
changed  only  in  that  p and  V become  pe  and  Ve  which  are  the  values  at  some  initial  condition 
where  the  stability  axes  are  fixed.  Lengths  are  divided  by  /o,  velocities  by  Ve,  time  derivatives 
by  Ve/lo  (for  each  differentiation)  and  mass  by  yeSlo  which  is  the  mass  of  a reference  volume. 
Aero-Normal  mass  is  called  p.  but  other  quantities  are  denoted  by  a superscript  bar. 


2.  CONSTANT  VELOCITY  UNDAMPED  SYSTEM 

The  equations  of  motion  of  a single  degree  of  freedom  system  may  be  expressed  in  terms 
of  Newton’s  equation  for  a force  and  a moment: 

force  = mass  x acceleration 

moment  = moment  of  inertia  x angular  acceleration. 

Ignoring  velocity  terms,  the  Z force  may  be  expressed 

Z = Za  a + Zo  (2.1) 

and  the  moment  as 

m = m„  a -f  mo.  (2.2) 

For  a body  of  mass  mt  and  moment  of  inertia  /»,  Newton's  equations  then  become: 

Z„  at  + Zo  = m,  1 

= 6. 


and 
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(2.3) 

(2.4) 


If  the  initial  condition  is  stable  and  primed  variables  are  changes  from  this  initial  condition, 
these  equations  may  be  written 

Za  a'  = m,  i (2.5) 

and  mx  a'  = Iy  6.  (2.6) 

Each  term  in  equations  2.5  and  2.6  can  be  made  non-dimensional  by  multiplying  and/or  dividing 
by  appropriate  constants,  as  follows: 

Za  , nte  lo 


IpeSlo  Ve2  Z 


which  is  written1**1 


Cta  «'=/!?. 


Similarly  for  the  moment  expression 


which  becomes 


mx  , Iy  lo*9 
lp.Vt*SI0  * “ kpeSlo3  V? 

Cm.  «'  = Iy  l 


This  non-dimensional  or  aero-normal  notation  is  used  throughout  this  work  except  where 
particular  cases  note  otherwise. 

These  equations  of  motion  (2.8  and  2.10)  are  connected  by  the  geometric  relation 

0 = « + i/K,  = « + i.  (2.11) 

This  leaves  two  independent  variables  with  two  separate  additive  solutions.  One  way  to  arrive 
at  each  solution  is  to  choose  appropriate  initial  conditions. 

2.1  Heave  Mode 

If  initial  conditions 

0 = 0 (2.12) 

* = 0 = 0 (2.13) 

and  i = K (2.14) 


are  imposed  on  the  system  of  equations  2.8,  2.10  and  2.11,  a trial  solution  can 

z = zo  + zi 
a = ao  + “l  e*,*'<+e, 

e = 00  + 01 

Equation  2.14  means,  in  conjunction  with  2.15,  that 

K = jo>z\  e}"t0)  = jotzi. 

Similarly  2.17  and  2.12  yield 

0 = 0o  + 0i 

and  2.17  and  2.13  give 

0 = ju>6\ 

2.19  and  2.20  further  imply  0O  = 0i  = 0 

so  0 remains  at  its  initial  zero  value.  Then  equation  2.11  becomes 

« = —i 

and  hence 

00  + 01  = -( K/V. y-4 

which  means,  after  evaluating  this  expression  for  t not  zero  and  t = 0,  that 

o!  = -K/V,  = -R 

and  ao  = — «i  e*. 


be  tried : 
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FIG.  1.  FORCE  AND  MOMENT  AXES  SYSTEM 
(Sine  definitions  of  angles) 


In  equation  2.10,  the  right  hand  side  term  0 is  zero  but  the  left  hand  side  contains  no  zero 
terms.  The  “heave”  mode  must  then  assume  infinite  moment  of  inertia  Iv , and  for  a constant 
velocity  system  equation  2.8  becomes: 

K u Iq 

- C*„  — e>"-‘  = ^ jwK  el“‘.  (2.25) 

Cancelling  terms, 

jw  = — (V,  Cix)/(n  l0),  which  is  real,  = A,  say.  (2.26) 

The  heave  mode  is  now  seen  to  be  an  exponential  response  with 

z = zo  + (AT/A)  e"  (2.27) 

<*  = -K  = -(K/Ve)  e“  (2.28) 

where  A = — Veftn  /„)  Cla. 

2.2  Pitch  Mode 

Taking  new  initial  conditions 

z(0)  = 0 
9(0)  = K 

and  0(0)  = a(0)  = 0 

then  equations  2.15  to  2.17,  the  same  trial  solution  as  before,  give 


(a)  z(0)  = 0 = jwzi  e>“°  (2.30) 

z = z0  + zi  cos  u>t  (2.31) 

and  a = ai  sin  wt  (2.32) 

and  6 = 6i  sin  wt.  (2.33) 

(b)  6 = w8 1 COS  wt 

8(0)  = w8t  = K (2.34) 

0i  - K/w 

8 = (K/w)  sin  wt.  (2.35) 

Substituting  back  in  our  equations  of  motion, 

2.10  -+  Cmx  «i  sin  wt  = — Ty  (K/w)  (/o2/k«z)  w 2 sin  wt 

Cma  ai  = -Ty(K/w)w*.  (2.36) 

In  equation  2.8,  the  right  hand  side  term  2 is  zero  but  the  left  hand  side  contains  no  zero 
terms.  The  “pitch”  mode  must  then  assume  infinite  aero-normal  mass  /x.  This  further  implies 
that 

ai  = 0i  = K/w.  (2.37) 

Substituting  for  <*i  in  2.36  and  clearing, 

<=2  = -CmJlv.  (2.38) 


The  pitch  mode  is  thus  a steady  sinuosoidal  oscillation  in  both  pitch  and  incidence. 

The  system  of  this  section  is  not  physically  realizable.  Real  systems  have  losses  (damping) 
and  in  the  next  section  damping  terms  will  be  considered. 

3.  THE  DAMPED  CONSTANT  VELOCITY  SYSTEM 

The  equations  of  motion  considering  velocity  terms  in  the  aerodynamic  forces  become 

Cta  a -f  Ctjl  S'  + Czq  q — (3-1) 

and  Cm<1 «’  + Cm^  * + Cm,  q = ly  tf  (3.2) 
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with 


0 = a'  + z 


(3.3) 


where  primed  variables  are  changes  from  initial  conditions  and  all  terms  are  aero-normal. 
The  same  two  sets  of  initial  conditions  are  considered  for  this  system  as  for  the  undamped 
system  in  the  previous  section. 


3.1  Heave  Mode 

Assume  initial  conditions 

0 = 0 i 

8 = 0 J>  (3.4) 

As  in  section  2.1,  the  first  two  conditions  are  met  by 


and  a solution  of  the  form 

0(0  = 0 

(3.5) 

z = zo  + zi  eH 

(3.6) 

is  tried. 

a = ao  + ai  eA< 

(3.7) 

Equation  3.3  gives 

0 = ai  eA‘  + (zi  AeA‘)/F« 

(3.8) 

at  = — Zl  S. 

(3.9) 

From  3.4, 

z(0)  = zi  A eA0  = K 

= zi  A 

.'.  zi  = K/X  or  zi  = R/X. 

(3-10) 

Hence  from  3.9  and  3.10, 

ai  = —K. 

Substituting  back  in  equation  3.1, 

(3-11) 

—K  CZx  eA(  - KJ  CZi  <?A<  = n (K  X)  eA( 

(3.12) 

so  A = - CzJ(p  + Cz.).  (3.13) 

3.2  Pitch  Mode 

Assuming  new  initial  conditions 

z(0)  = 0 

9(0)  = K 

> (3.14) 

and  trial  solutions 

0(0)  = a'(0)  = 0 

z = zt  ] 

a'  = 0 = ao  eA<  e,utl  ; 

l (3.15) 

then 

q(0)  = K = ao  AeA0  e)a,°  + ju>  ao  eA0  ei<j0  (3.16) 

and 

0(0)  = 0 = ao  eA0  e*ul 

= ao  sin  (<u0) 

(3.17) 

as  the  cosine  term  must  be 

zero  also. 

From  3.17 

a'  = 0 = ao  eA<  sin  (<uf ) 

(3.18) 

and  3.16  becomes  K = 

ao  A 0 -f-  aof°a)  COS  (a>0) 

= 

tu  ao 

ao  = 

K/w 

(3.19) 
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ami  the  equation  of  motion  becomes  3.2 

Cm ,a  (K/w)  eM  sin  (cut)  + (Cmj  + C,n<|)  {[(A3)/<y]  e A(  sin  (cut)  + K ext  cos  (cot)}  = 


= Iy  K {23  cos  (to/)  + [(32  — ffi2)/dj]  sin  to/}  ext.  (3.20) 

By  equating  the  sine  and  cosine  terms,  and  clearing, 

Cinj  + Cm,  = 23  Iy  (3.21) 

and  Cm,  “F  (Cm±  “l-  CmQ)  ^ — Iy  (32  — o2).  (3.22) 

Substituting  3.21  into  3.22  then  gives 

Cm,  = - (32  + *2)  h (3.23) 

= — ao2  Iy-  (3.24) 


This  is  of  the  same  form  as  in  the  previous  section  (Eqn  2.38)  but  to  is  now  the  complex 
quantity  too.  For  a dynamically  (as  well  as  statically)  stable  system,  the  real  part  of  <o0  must  be 
negative  If  positive,  i.e.  A > 0,  any  disturbance  in  pitch  will  grow. 

The  two  modes,  heave  and  pitch,  of  3.1  and  3.2,  exist  concurrently,  and  a suitable  choice 
of  initial  conditions  shows  both.  Flowever,  as  long  as  the  time  scales  of  the  two  disturbances 
are  sufficiently  different,  say  one  more  than  forty  percent  longer  than  the  other,  the  resultant 
remains  simply  the  sum  of  the  two  disturbances. 

The  damped  pitch  disturbance  occurs  in  aircraft  (fixed  controls)  and  can  be  simulated  in 
a wind-tunnel  with  low  friction  pivots.  The  heave  mode  does  not  however  occur  in  this  constant 
velocity  form.  In  flight,  speed  will  change  as  the  aircraft  climbs  or  falls  modifying  the  result. 
In  a wind  tunnel  it  is  difficult  to  simulate  translations;  a rotation  about  a point  not  at  the  reference 
centre  will  be  considered  instead. 

4.  SPRING  RESTRAINED  PITCH  SYSTEM 

Consider  a system  where,  as  well  as  the  aerodynamic  terms,  there  are  structural  stiffness 
(Kg)  and  damping  (Z^)  terms  but  heaving  motions  are  constrained.  For  the  wind-off  or  tare  case. 


but  non-dimensionalizing  as  for  the  wind-on  case  to  be  considered  next, 

IyS  ^ K0  d + Kg  S (4.1) 

which  has  the  solution 

6 = 60  (4.2) 

with  Kg  = 2 3(  Iy  (4.3) 

and  Kg  = — (cut2  + \2)  /„  (4.4) 

where  3<  = A,  l0/  V,  and  = wt  lol  Ve 

with  V,  referring  to  the  wind-on  case. 

Since  a pitch  system  with  a = 6 is  being  considered,  the  wind-on  case  becomes; 

Iy  8 — (Kg  + Cm,)  ft  + (Kg  + Cnij  + Cm,)  8 (4.5) 

which  has  the  solution 

0=0oe<  A +/<*>(  (4.6) 

with  Kg  + Cmi  + Cm,  = 23  Iy  (4.7) 

and  Kg  + Cma  = - (a>2  + X2)  Iy.  (4.8) 

From  equations  4.3  and  4.7 

Crai  + Cm,  = 2(3  - I,)  Iy  (4.9) 

and  from  equations  4.4  and  4.8 

Cm,  = — (A2  — 3(2  + a2  — at2)  Iv  (4.10) 

= (®n,2  — g>o2)  Iy  (4.11) 

using  the  results  of  Section  3. 
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j. 


The  terms  on  the  left  in  4.9  and  4. 10/1 1 are  the  same  as  those  responsible  for  the  behaviour 
of  the  pitch  system  in  Section  3 and  the  equations  are  the  difference  between  two  sets  of  results 
of  the  form  of  that  Section.  This  convenient  relationship  means  that  this  system  is  often  used 
in  wind-tunnel  testing. 

5.  SPRING  RESTRAINED  PITCHING  AND  PLUNGING  SYSTEM 

If  the  centre  of  rotation  is  moved  from  the  origin  as  in  Section  4 to  (—  xo,0, 0)  in  the  body 
co-ordinates,  then 

6 = * + t (5.1) 

= « + HVt 

and  9 = z/x0.  (5.2) 

The  equation  of  motion  is  far  more  complicated 
PXO'  A - ly  ^ ~ Rg  B -f-  Rg  § f ( ni,  a -f  Cf 71^  & -f-  Cm,  I Cj„  .Xo  oc  Cz^  Xo  oc  • Cj,  Xo  (5.3) 
Equations  5.1  and  5.2  give 

B = x + x0§.  (5.4) 

The  moment  of  inertia  about  the  new  yi  axis 

A/,  = A/ 4- Mo2.  (5.5) 

Substituting  5.4  and  5.5  in  equation  5.3  g'^es 

Ai,  & — (Re  4"  Cma  4*  xo  Cz3 ) 8 -|-  (Rg  xo  f 4-  Cm,  4-  Cm,  — xo2  CZx  4-  xo  CZi  + xo  CZq)  5 — 

- xo  (Cmi  4-  xo  Czi)  If.  (5.6) 

Tare  values  give 

Re  = 23,  A,,  (5.7) 

and  Re  = - (a>,2  4-  V)  A/,  = — ®o,2  A/,-  (5.8) 

For  the  wind-on  case,  try  a solution 

0 = 0o<?<a+M  (5.9) 

Then  Rg  -f-  Cm,  4-  Cm,  • xo  (C^  4-  C»,  Cm^)  — xo2  CZl  = 2 3 (7Vl  4-  xo  Cm,  4-  xo2  C^)  (5. 10) 

and  Rg  4-  C*a  + x0  CZx  = — (32  + a2)  (/y,  4-  x0  Cmi  + x02  Cfi)  = 

= — coo2  (A i,  4-  xo  Cmj  4"  xo2  CZjt).  (5.1 1) 

Substituting  5.8  into  5.1 1, 

Cm,  4“  Xo  Ct3  = — (too2  — ®o,2)  Ac,  — ©o2  xo  (Cmi  4-  xo  CZi).  (5.12) 
Similarly  5.7  into  5.10  gives 

(Cm,  4*  xo  CZi)  4 - (Cm,  •-  xo  Cz,)  — xo  (Cmj  -f-  xo  Cj,)  = 2 (X  — 3,)  JVi  + 23  xo  (Cmi  4-  xo  C2j). 

(5.13) 

Using  5.12,  equation  5.13  becomes: 

(Cfflj  4-  xo  C*,)(l  4 xo2  coo2  — 23 xo)  4-  (Cm,  4-  xo  Cz,)  = [2(3  — 3,)  — xo  (&o2  — coo,2)]  /y,.(5.14) 

This  is  complicated.  However  if  values  for  xo  = 0 are  known  it  allows  measurement  of 
Cz<1  and  estimation  of  the  split  between  Cra,  and  Cm,  which  have  appeared  together  up  to  now. 

6.  FORCED  SPRING-RESTRAINED  PITCH  SYSTEM 


To  the  system  of  Section  4 a torque  7\>  e,ut  is  added.  The  tare  case  becomes: 


1 y $ — fifi  0 Kf)  0 

i To  eW. 

(6.1) 

The  homogeneous  solution  (To  0)  was  given  in 

0 = 0O 

Section  4 

(4.2) 
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where  Re  = 2A,  lv  (4.3) 

and  Re  — — (®(2  + ^i2)  Jy  (4.4) 

= - ®0,2  ly.  (6.2) 

The  complete  solution  in  the  tare  case  is  thus 

0 = 0 o -f  6i  (6.3) 

where  it  may  be  shown  that 

£ = tan-1  [(—Rg  &)/ (ly  ft2  + Re)],  which  using  4.3  and  4.4, 

L-  = tan"1  [(-2A,  ft)/(ft2  - so,2)]  (6.4) 

T 

and  also  6l  = - + (IP  - so,2)2]*'  (6,5) 


Note  that  for  zero  tare  damping,  i.e.  At  = 0,  then  £ = 0 so  that  the  oscillation  is  in  phase  with 
i he  applied  torque,  and 


0i  = To/Uy  (S0,2  - ft2)] 

(6.7) 

which  is  infinite  at  resonance. 

The  wind-on  equation  of  motion  is 

IyB  = 

Re  0 T"  Rd  $ + Cmx  0 + Cmi  ^ Cm,  S + To  elot. 

(6.8) 

The  solution!  is  of  the  form 

6 = 0O  ea+iu)t  + 01  et<fl<+£) 

(6.9) 

where 

Cma  = - (so2  - ®o,2)  Jv 

(6.10) 

Cmi  Cm,  — 2(A  — A,)  Iy 

(6.11) 

and 

£ = tan-1  [(2A  ft)/(s02  - ft2)] 

(6-12) 

remembering 

So2  ==  S2  ~h  A2 

(6.13) 

and  where 

01  = To/(/y  [4A2  ft2  + (ft2  - ®o2m- 

(6.14) 

For  zero  damping  (A  = 0)  equation  6.14  reduces  to 

01  = To/[/y  (mo2  - ft2)] 

(6.15) 

which  again  gives  infinite  amplitude  for  the  resonance  condition  a»o  = ft. 

For  non-zero  damping,  the  resonance  condition  gives  £ = tan-1  oo  = 

90  degrees  for  both 

tare  and  wind-on  cases.  In  the  tare  case, 

ft t = So, 

(6.16) 

then 

01  = To/(2/y  A,  ft). 

(6.17) 

For  wind-on 

ft  = So 

(6.18) 

then 

01  = To/(2lv  A ft). 

(6.19) 

Combining  these  equations 

Cmi  + C„,  = T0/(O  6l)  - To, /(ft,  0i,).  (6.20) 

If  tare  and  wind-on  tests  have  the  same  amplitude,  and  stiffness  can  be  varied  to  obtain  the  same 
frequency,  then,  providing  Ke  is  the  same, 

Cnk  + Cm,  = ( To  - fo,)/(fl  0l).  (6.21) 

This  is  similar  to  the  results  obtained  for  free  oscillation  with  the  decay  rate  replaced  by 
a torque  input.  It  has  an  advantage  in  that  the  output  is  steady  and  measurements  can  be  averaged 
over  long  time  periods  for  high  accuracy.  As  with  the  free  system  an  estimation  of  the  split  of 
Cmi  and  Cn,  requires  multiple  axis  locations. 


t S.  Timoshenko,  Vibration  Problems  in  Engineering.  Van  Nostrand  1928. 
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Test  section 


(a)  Relations,  angles  and  lengths  for  towed  models  (h^=— x 6 + I 0C) 


Incidence 


(b)  Relations,  angles  and  velocities  for  towed  models  (a  = 6 — z/v) 


(c)  Forces  and  moment  Towed  model  (T  = — x cos  (6  + 6C)  — Z sin  ( 6 + 6C ) 

+ mx  sin  (0  + 01  — W sin  0. 


FIG.  2.  RELATIONSHIPS  REQUIRED  FOR  SECTION  7 
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I 7.  CABLE  TOWED  MODELS 

7.1  Pitch/Incidence  System 

| Consider  a wind  tunnel  model  on  a long,  flexible  cable.  The  constraint  conditions  are  taken 

to  be: 

(i)  The  cable  is  inextensible  (length  /),  massless  and  dragless; 

(ii)  No  moment  about  the  point  of  attachment  is  provided  by  the  cable  which  is  attached 
distance  x ahead  of  the  model  centre  of  gravity; 

(iii)  No  disturbance  in  roll  or  yaw  exists  implying  that  the  roll-yaw  system  must  be  statically 

> and  dynamically  stable  and  have  natural  frequencies  sufficiently  different  from  those  of 

^ the  pitch/incidence  system  to  minimize  cross-excitation  of  one  system  on  the  other. 

Condition  (i)  means  that  the  force  on  the  cable  must  equal  the  tension  T in  the  cable, 
i Resolution  of  all  the  forces  acting  leads  (Fig.  2a)  to  the  expression 

T = — A' cos  ( 8 + 8C)  - Z sin  (0  + 8C)  + mx  sin  ( 8 + 8e)  - sin  8e  (7.1) 

for  the  tension  in  the  cable,  and 

—A L = Z cos  8 + IF  — X sin  8 -(-  7"  sin  8e  (7.2) 

for  the  force  unbalance  in  the  vertical  direction.  This  will  equal  zero  for  a static  equilibrium. 
Similar  resolution  yields 

Aw  = m +Tx  sin  ( 6 + 8C)  (7.3) 

for  the  resulting  pitching  moment.  In  each  of  these  three  equations,  the  terms  have  not  been 
made  non-dimensional  as  yet,  in  order  to  show  the  precise  nature  of  some  approximations  made. 
For  8 and  8C  both  small,  from  Fig.  2, 

T sin  8C  » [-x  -Z  (8  + 6C)  + mx  (8  + 8C)  - W 8C]  8C 


* -X8C  (7.4) 

being  first  order  in  8C,  and  similarly 

Tx  sin  (0  + 8C)  * -Xx  (8  + 8C).  (7.5) 

The  nett  vertical  force  is  then 

-A L x Z + W - X(8  + 8C)  (7.6) 

and  the  nett  moment  about  the  y-axis  is 

Aw  « w — Xx  (8  -f  8C).  (7.7) 

Now  8,  8C,  z and  a are  related  (Fig.  2b)  by 

z = x8-l8c  (7.8) 

and  a.  = 8 — z/V.  (7.9) 


From  equations  7.8  and  7.9  it  is  clear  that  there  are  still  only  two  independent  variables  as 
in  Sections  2 and  3.  So  comparisons  can  be  made  with  the  equations  derived  in  those  sections, 
the  equations  of  motion  are  considered  in  non-dimensional  form. 

If  the  initial  conditions  represent  equilibrium,  then 


Ctq  + W = 0 (“lift”  = weight) 

(7.10) 

and 

Cm0  + Cm? . 7)  = 0 (zero  moments) 

(7.11) 

and  the  equations  of  motion  become 

CzI  a + Cu  4 + C2,  q — (8  + 8C)  (Cr0  + Czx  « + CXi  4 + CXq  q)  — ft  z 

(7.12) 

and 

Cj»j  « + Cmj  4 + Cm,  q — x(8  -f  8e)  (C.r„  + CXx  a + CXi  i + CZq  tj)  = Iy  9. 

(7.13) 

It  is  assumed  that  the  solution  of  these  equations  can  be  divided  into  two  additive  modes 
as  in  Sections  2 and  3.  It  is  further  assumed  that  one  way  of  representing  these  modes  is  by 
separating  Equation  7.9  by  zeroing  one,  then  the  other,  of  the  right  hand  terms.  These  modes 
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are  then  the  “pitch”  and  "heave"  modes  of  Sections  2 and  3,  and  the  equations  of  motion 


become: 

- Cz,  f - Cz,  ! - (f/7)  ( -Cx0  + Ct,  5 + Cz,  *)  = fi-  (7. 14) 

for  the  mode  with  z = —a.  (7.15) 

and  C*,  0 + (Cm,  + Cm,)  8 - X0  [Cz,  + CXJ  + ( Cz , + Cx,)  <*]  = /*  £ (7. 16) 

for  the  mode  with  0 = a.  (7.17) 

Both  7.14  and  7.16  are  non-linear  differential  equations.  The  effect  of  the  tow  cable  is  to 
introduce  a non-linear  additional  stiffness  to  each  mode.  If  terms  of  more  than  first  order  can  be 
ignored,  these  equations  become: 

- Cz,  z - CXi  i + Cx>  ( z/1 ) = M * (7. 1 8) 

and  (Cm,  - * Cz,)  6 + (Cm,  + Cm,)  8 = /„  8.  (7.19) 

Equation  7.19  is  identical  in  form  to  the  spring  restrained  pitch  system  of  Section  4 with 
spring  constant 

Ke  = -x  Cz.  (7.20) 

and  spring  damping 

Ke  = C.  (7.21) 


8.  FORCES  AND  MOMENTS  ON  A RIGIDLY-DRIVEN  MODEL 

If  a model  of  mass  me  and  of  moment  of  inertia  about  its  centre  of  mass  /„  is  driven  rigidly, 
then  a force  and  moment  balance  between  the  model  and  the  driving  mechanism  measures 
inertial  as  well  as  aerodynamic  terms.  The  balance  axis  system  is  a body-axis  set,  and  not  the 
fixed  “stability  axes”  used  for  other  methods  of  driving  oscillatory  models.  Two  positions  of 
the  axis  of  rotation  are  considered  as  well  as  a simple  translation.  The  balance  is  considered 
to  be  rigid  in  comparison  with  the  movement  applied. 


8.1  Centre  of  Rotation  at  Reference  Centre  (0,  0,  0) 

For  a rigid  sinusoidal  pitch  oscillation  about  the  y-axis, 

a = ao  + sin  wl  (8.1) 

6 = do  + 8i  sin  cut  (8.2) 

a = 8 — a 0i  cos  a>t  (8.3) 

U = 8 = — o2  0i  sin  <ot.  (8.4) 

The  z force  coefficient  has  the  form,  where  g is  the  acceleration  due  to  gravity  and  /x  is  the 
non-dimensional  model  mass  (m«), 

Cz  = (i  g cos  (ao  + 0i  sin  cot)  -f  Cz,  + Cz,  0i  sin  cot  -H  ( Ct i + Clq)  ® 0i  cos  cot.  (8.5) 
For  a small  amplitude  0i  this  may  be  approximated  by 


Cz  — (jig  cos  ao  + CzJ  + (Cz,  — fig  sin  ao)  0i  sin  cot  + (Cz,  + C2.)  to  0i  cos  cot  (8.6) 

where  the  three  terms  are  the  static,  in-phase  and  quadrature  components  of  the  force  with 
respect  to  the  motion,  respectively. 

The  pitching  moment  coefficient  is  given  by 

Cm  = (Cm,  4-  Cm,i)  4-  (Cm,  - O*  /„)  #1  sin  cot  + (Cm,  + Cm,)  O »l  COS  cot  (8.7) 

where  the  same  comment  about  the  three  terms  applies  and  where  »/'  is  a fixed  control  deflection. 


8.2  Centre  of  Rotation  Shifted  Along  the  x-axis  (xo,  0,  0) 

If  the  model  is  rotated  about  an  axis  xo  ahead  of  the  y-axis,  a and  0 are  no  longer  co-incident 
quantities.  As  always,  from  the  definitions  of  the  angles  of  pitch  and  incidence. 


In  this  case 


0 = * 4-  i,  i.e.  = a 4-  z/Vt. 
0 = r/x o,  also. 


(8.8) 

(8.9) 
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Thus  0 = a + X oq  (8,10) 

which  is  the  same  as  equation  5.4  for  the  spring-restrained  system. 

The  model  now  has  moment  of  inertia 

Iyi  = Iy  + X02  me 

IVl  = Iy  + X„2  (8.11) 

The  oscillating  model  now  has  a linear  acceleration  term  in  its  force  equation.  The  equations 
of  motion  are  now 

0 = 0o  +0i  sin  <ot  (8.12) 

and  a = 0O  + 0j  sin  cut  — xo  a>  0i  cos  tut  (8.13) 

with  derivatives 

8 = q = 0i  © cos  u>r  (8.14) 

a = 0i  o cos  tut  + xo  0i  a2  sin  tot  (8.15) 

1 = — xo  ©2  0i  sin  tut  (8.16) 

and  8 = q = — 0i  u>2  sin  <ut.  (8.17) 

Then  the  z force  coefficient  has  the  form 

Ct  = C*0  + fi  g cos  0 + ii  z + Clx  a + Cti  *■+  C,t  q.  (8.18) 

If  the  oscillation  amplitude  0i  is  small 

cos  0 a cos  0o  — sin  (0o)  0i  sin  <ut.  (8.19) 

Including  all  the  constant  terms  in  Cz0  to  make  Cz>'  and  substituting  for  the  angles  in  8.18, 

Ci  = <V  + [Cia  + fx  (xo  a2  — g sin  0O)  + Cl.  xo  o2]  0i  sin  tut  + 

+ (Cii  + Ctt  — xo  CtJ  as  0i  cos  <ut.  (8.20) 

The  pitching  moment  is  even  more  complicated.  The  moment  about  the  new  axis  is  given  by 

m'  = m + xoZ  (8.21) 

i.e.  Cm>  = Cm  xo  Ci.  (8.22) 

The  balance  however  measures  moments  about  the  reference  axis,  i.e.  Cm.  The  expression  *or 
this  is  exactly  as  before: 

Cm  — (Cm„  + Cm<i  t)')  + (Cma  — ®2  fy)  61  sin  tut  + (Cmi  + Cm,)  ® 01  COS  tut.  (8.7) 
Moments  about  the  new  axis  can  be  obtained  using  Equations  8.7,  8.20  and  8.22. 

8.3  Rigid  Translation  Along  the  z-axis 


If  the  motion  inexorably  impressed  on  the  model  is  a z translation  and  not  a rotation, 


pitch  angle  0 remains  a constant,  say  zero.  Then  from  equation  8.8 

« = -i  (8.23) 

and  for  a motion  z = zo  + zi  sin  tut  (8.24) 

a = — zi  a cos  tut  (8.25) 

and  i ==  zita2  sin  tut.  (8.26) 

The  z-force  and  pitching  moment  coefficients  are  given  by,  ignoring  static  valves. 

Cl  = Cia  « + Ctit  + M I (8.27) 

a*>d  Cm  — Cm,  « + Cmj,  £ (8.28) 

and  using  8.25,  8.26  these  become 

Ci—  — Cia  zi  tu  cos  tut  + (Cid  — fi)  zi  a2  sin  tut  (8.29) 
Cm  = — Cm,  zi  as  cos  tut  + Cmi  zi  o*  sin  tut.  (8.30) 
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and 


9.  CONCLUSION 


The  equations  of  motion  have  been  developed  for  a simplified  free  flight  pitch-yaw  system 
and  for  spring-mounted  and  rigidly-driven  wind  tunnel  systems.  For  the  spring-mounted  system 
both  initial  displacement  and  forced-oscillation  conditions  have  been  examinee.  A simplified 
cable  towed  system  has  also  been  derived.  These  equations  will  be  used  in  later  Notes  of  this 
Series. 
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